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Abstract 

We present the spectra of negative charged particle production in Au+Au collisions 
from STAR for the hrst phase of the RHIC Beam Energy Scan Program measured over a 
wide range of collision energy 7-200 GeV, and transverse momentum of produced 

particle in different centralities at \ri\ < 0.5. The spectra demonstrate strong dependence 
on collision energy which enhances with pt- An indication of self-similarity of negative 
charged particle production in Au+Au collisions is found. The constituent energy loss as 
a function of energy and centrality of collisions and transverse momentum of inclusive 
particle was estimated in the 2 ;-scahng approach. The energy dependence of the model 
parameters - the fractal and fragmentation dimensions and "specihc heat was studied. 
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I. INTRODUCTION 


A major goal of high-energy heavy-ion collision experiments at RHIC is to determine 
the phase diagram for nuclear matter and describe its features in framework of 
the Quantum Chromodynamics (QCD). The most experimentally accessible way to 
characterize the QCD phase diagram is in the plane of temperature (T) and the baryon 
chemical potential (hb) [H- 

Heavy-ion collisions at RHIC have provided evidence that a new state of nuclear matter 
exists at the top RHIC energy ^/snn =200 GeV |2_ 
a suppression of particle production at high [6|-l9|, a large amount o 
number-of-constituent-quark (NCQ) scaling of V 2 at intermediate Pt 


This new state is characterized by 
liptic flow (^ 2 ), 


correlated yields at large Ar/ and A0 ~ 0 (ridge) llMl3|. 


10| and enhanced 


Figure [H shows a schematic layout of the phases, along with hypothesized indications 
of the regions crossed in the early stages of nuclear collisions at various beam energies. 



Phc. 1: Phase diagram of nuclear matter 



The idea of the Beam Energy Scan (BES) Program at RHIC was to vary the 
collision energy with aim to search for the QCD critical point and QCD phase boundary 
The first phase of the Beam Energy Scan program (BES-I) along with the existing 


top energy at RHIC has allowed the access to a region of the QCD phase diagram (Fig. [T]) 
covering a range of baryon chemical potential {ps) from 20 to 420 MeV corresponding to 
And-An collision energies from ^^snn =200 to 7.7 GeV, respectively. Results from BES-I 
have further supported the evidence for the quark-gluon plasma (QGP) discovery at the 
top RHIG energy y/s^N =200 GeV [2^. The results of the search for the critical 
point and the first-order phase boundary have narrowed the region of interest to collision 
energies below y/sj^ =20 GeV. 

The principal challenge remains localization of a critical point on the QGD phase 
diagram. Near the QGD critical point, several thermodynamic properties of the system 























such as the heat capacity, compressibility, correlation length are expected to diverge with 
a power-law behavior in the variable e = {T — Tc) /T^, where is the critical temperature. 
The rate of the divergence can be described by a set of critical exponents. The critical 
exponents are universal in the sense that they depend only on degrees of freedom in 
the theory and their symmetry, but not on other details of the interactions. This scaling 


postulate is the central concept of the theory of critical phenomena 2ll,l22j. Current lattice 


QCD calculations suggest that key features of the phase diagram like the critical point 
and the hrst-order phase transition lie within the /is reach of the RHIC BES Phase-II 
program (see and refernces therein). 


II. EXPERIMENT AND DATA ANALYSIS 

The preliminary results presented here are based on data taken at STAR in Au+Au 
collisions at ^/s]^ = 7.7, 11.5, 39, 62.4 and 200 GeV in the year 2010 while 19.6 and 
27 GeV were collected in the year 2011. The transverse momentum spectra of negatively 
charged particle production in the mid-rapidity over a wide range of centrality are reported 
in this paper. 


STAR detector 


The main detector used to obtain the results on particle spectra is the Time Projection 
Ghamber (TPG) 23|. The TPG is the primary tracking device at STAR and can track up 
to 4000 charged particles per event. It is 4.2 m long and 4 m in diameter. Its acceptance 
covers ±1. units of pseudo-rapidity rj and the full azimuthal angle. The sensitive volume 
of the TPG contains PIO gas (10% methane, 90% argon) regulated at 2 mbar above 
atmospheric pressure. The TPG data are used to determine particle trajectories, momenta, 
and particle-type through ionization energy loss (dE/dx). STAR’s solenoidal magnetic held 
used for this low energy Au+Au test run was 0.5T. The Time of Flight (ToF) detector 


24| (with 27r azimuthal coverage and |? 7 | < 1.0) further enhances the PID capability. All 


events were taken with a minimum bias trigger. The trigger detectors used in these data 
are the Beam-Beam-Gounter (BBG) and Vertex Position Detector (VPD)[3- The BBGs 
are scintillator annuli mounted around the beam pipe beyond the east and west pole-tips 
of the STAR magnet at about 375 cm from the center of the nominal interaction region 
(IR), and they have an t] coverage of 3.8 < |? 7 | < 5.2 and a full azimuthal (27r) coverage. 
The VPDs are based on the conventional technology of plastic scintillator read-out by 
photomultiplier tubes. They consist of two identical detector setups very close to the 
beam pipe, one on each side at a distance of |14| = 5.6 m from the center of the IR. The 
details about the design and the other characteristics of the STAR detector can be found 
in 01 • 

All events selected for the analysis satished the cuts described below. The primary 








Table 1. Track selection criteria at all energies. 


|r/| DCA nFitPts nFitPts/nFitPoss ndE/dx pTpr/pTgl 
< 0.5 < 1 cm > 15 > 0.52 > 10 7/10 < & < 10/7 


vertex for each minimum bias event is determined by finding the best point of common 
origin of the tracks measured in the TPC. In order to reject events which involve 
interactions with the beam pipe and beam-gas interactions, the vertex positions along 
the X, y and longitudinal {z) beam directions are required to be less than 2 cm and 30 
cm, respectively. 

Track selection criteria for all data sets are presented in Table III Al Primary tracks 
are placed on the distance of closest approach (DCA) between each track and the event 
vertex. Tracks must have at least 15 points (nFitPts) used in track fitting. To prevent 
multiple counting of split tracks, at least 52% of the total possible fit point are required 
(nFitPts/nPitPoss). A condition is placed on the number of dE/dx points used to derive 
dE/dx values (ndE/dx). The cut for the pr ratio of primary and global tracks (pTpr/pTgl) 
were used to delete fail tracks. The results presented here are for transverse momentum 
Pt > 0.2 GeV/c and pseudo-rapidity \p\ < 0.5 were selected for the analysis. 

The efficiency is determined by measuring the reconstruction efficiency of Monte Carlo 
tracks embedded in real events and then run through the STAR reconstruction chain 
including realistic STAR geometry. For a given collision energy, centrality region, and 
particle species the efficiency is fitted as a function of pt using the functional form eff = 
A ■ exp{ — {B/pT)^}- The parameters A, R, and C are determined for each centrality for 
unidentified particles using efficiencies for identified particles. 

Figure 2 demonstrates dependence of reconstruction efficiency on transverse 
momentum of particle at energy ^Jsnn = 7.7 and 200 GeV and various centralities. 



a) b) 

Phc. 2: Reconstruction efficiency of negatively charged particles produced in Au+Au collisions 
at energies ^/snn = 7.7 and 200 GeV in the range \r]\ < 0.5 and various centralities as a function 
of Pt- 



















Centrality classes in Au+Au collisions at ^^snn = 7.7 — 200 GeV are defined using the 
number of charged particle tracks reconstructed in the main TPC over the full azimuth 
and pseudorapidity \ri\ < 0.5. This is generally called "reference multiplicity (refmult)"in 
STAR HQ- For each energy, a refmult "correction"is applied on the standard definition. 
This is based on the following steps: identifying bad runs, applying acceptance/efficiency 
correction of refmult for different vertex positions, and performing low refmult correction 
for trigger inefficiencies for different 2 ;-vertices. The refmult distribution is compared/fitted 
to Monte-Carlo Glauber model simulations. The centrality is defined by calculating the 
fraction of the total cross-section obtained from the simulated multiplicity. The events are 
divided into following centrality classes: 0 — 5%, 5 — 10%, 10 — 20%, 20 — 40%, 40 — 60%, 
60 — 80%. Table III Al shows the total number of events that are used for the analysis for 
each energy after the above event selection cuts. 


Table 2. Total events analyzed after all analysis cuts for various energies. 


Year 

(GeV) 

Events 

(Millions) 

7.7 

2010 

2 

11.5 

2010 

7 

19.6 

2011 

17 

27 

2011 

33 

39 

2010 

108 

62.4 

2010 

62 

200 

2010 

35 


III. RESULTS AND DISCUSSION 

A. Transverse momentum spectra 

The transverse momentum spectrum of hadrons produced in high-energy collisions 
of heavy ions reflects features of constituent interactions in the nuclear medium. The 
medium modification is one of the effects (recombination, coalescence, energy loss, multiple 
scattering,...) that affects the shape of the spectrum. The properties of the created medium 
are experimentally studied by variation of the event centrality and collision energy. 

Figures [3] and 0] shows the negatively charged hadron yields in Au+Au collisions at 
y/SNN =7.7, 11.5, 19.6, 27, 62.4 and 200 GeV at mid-rapidity \ri\ < 0.5 as a function 
of transverse momentum pt- The results are shown for the collision centrality classes of 
0 — 5%, 5 — 10%, 10 — 20%, 20 — 40%, 40 — 60%, 60 — 80%. The distributions are measured 
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Phc. 3: Transverse momentum spectra of negatively charged particles produced in Au+Au 
collisions at energies y^siviv =7.7, 11.5, 19.6 and 27 GeV in the range \r]\ < 0.5 and for different 
centralities (0-5%,5-10%,10-20%,20-40%,40-60%,60-80%). Errors shown are statistical only. 


over a wide momentum range 0.2 < < 12. GeV/c. The multiplied factor of 10 is used 

for visibility. As seen from Figs. [3] and 0] spectra fall more than eight orders of magnitude. 
Spectra demonstrate a strong dependence on collision energy and centrality at high px- 
The exponential behavior of spectra at low px and the power behavior of spectra at high px 
are observed. Figure 111(d) shows the transverse momentum spectra of negatively charged 
particles produced in Au+Au collisions at energies y/s^N = 7.7, 11.5, 19.6, 27, 39, 62.4 
and 200 GeV in the range \p\ < 0.5 for most central events (0 — 5%). The difference in 
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c) d) 

Phc. 4: Transverse momentum spectra of negatively charged particles produced in Au+Au 
collisions at energies =39 (a), 62.4 (b) and 200 GeV (c) in the range \r]\ < 0.5 and for 

different centralities (0-5%,5-10%,10-20%,20-40%,40-60%,60-80%). (d) Spectra for all energies 
y/s]\fi\f = 7.7 — 200 GeV in central events (0-5%). Errors shown are statistical only. 

particle yields between different energies increases with transverse momentum It is about 
three orders of magnitude for 0 — 5% Au+Au collisions at pt = 5 GeV/c between energies 
y/SNN =11.5 and 200 GeV. As we can see the data presentation in terms of dimensional 
quantities (a particle yield and transverse momentum) demonstrates a strong dependence 
on the scale. 











B. Self-similarity of negative particle production 


The study of critical phenomena in nuclear matter is most preferable in terms of 
dimensionless variables. The idea of the Beam Energy Scan programs at RHIC is to vary 
the collision energy and look for the signatures of QCD phase boundary and QCD critical 
point i.e. to span the phase diagram from the top RHIC energy to the lowest possible 
energy. Turn-off of signatures of QGP established at ■^/SNN =200 GeV would suggest the 
crossing of phase boundary. 

Self-similarity of hadron interactions is assumed to be a very fruitful concept to study 
collective phenomena in the hadron and nuclear matter. Important manifestation of this 
concept is a notion of scaling itself. Scaling in general means self-similarity at different 
scales. The approach based on 2 ;-scahng suggected in j^-30| is treated as manifestation 
of the self-similarity property of the structure of colliding objects (hadrons, nuclei), the 
interaction mechanism of their constituents, and the process of constituent fragmentation 
into real hadrons. We use the approach for analysis of the obtained spectra (see Figs. I3K|) . 
The approach relies on a hypothesis about self-similarity of hadron interactions at a 
constituent level. The assumption of self-similarity transforms to the requirement of 
simultaneous description of transverse momentum spectra corresponding to different 
collision energies, rapidities, and centralities by the same scaling function 'ip{z) depending 
on a similarity parameter z. The scaling function for the collision of nuclei is expressed 
in terms of the experimentally measured inclusive invariant cross section (Ed^a/dp^), the 
multiplicity density {dN/dp), and the total inelastic cross section (uin), as follows: 
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(dN/dp) (Tin 




( 1 ) 


Here J is the Jacobian for the transformation from {p^,y} to {z,p}. The function ip is 
interpreted as a probability density to produce an inclusive particle with the corresponding 
value of 2 ;. The similarity parameter z is expressed via momentum fractions (xi, X 2 , Pa, Vb), 
multiplicity density, and parameters 61 ^ 2 , ^a,b, c, as follows: 


z = ZqG 


( 2 ) 


where zq = ^/sT/[{dNch/dp\oYmN], G = (1 - xi)'^i(l - X 2 )‘^"(l - |/a)'“(l - VbY’’- 
Here is a nucleon mass. The quantity y/s± is the transverse kinetic energy of 

the binary subprocess consumed on the production of the inclusive particle with mass 
Till and its counterpart with mass m 2 . The constituents of the incoming nuclei carry 
fractions Xi,X 2 of their momenta. The inclusive particle carries the momentum fraction 
Ha of the scattered constituent. The parameters 61^2 and ea,b describe structure of the 
colliding nuclei and fragmentation process, respectively. The parameter c is interpreted as 
a "specihc heat "of the created medium. Simultaneous description of different spectra 
with the same puts strong constraints on the values of these parameters, and 

thus allows for their determination. It was found that 61^2 and c are constant at high 








Phc. 5: Scaling function 'ip{z) as a function of self-similarity parameter 2 for corresponding 
transverse momentum spectra of negatively charged particles produced in Au+Au collisions at 
energies y/s^N = 7.7 — 200 GeV in the range \r]\ < 0.5 for most central events (0-5%). 


energies, ea,b depends on mnltiplicity and mi = m 2 = m. For the obtained values 
of hi, 2 , ^a,b and c, the momentum fractions are determined to minimize the resolution 
Q~^{xi,X 2 ,ya,yb), which enters in the definition of the variable z. The system of the 
equations dfl/dxi = dQ/dx 2 = dfl/dya = dfl/dyb = 0 was numerically resolved 
under the constraint {xiPi -|- X 2 P 2 — p/yaY = [xiMi -|- X 2 M 2 + m/r/^ )^, which has sense 
of the momentum conservation law of a constituent subprocess 28l-l30|. The extended 
version for the approach was applied for description of hadron yields in nucleus-nucleus 
collisions j^. The fractal dimension of the nucleus is expressed in terms of the nucleon 
fractal dimension (h) and the atomic number of nucleus (A) as follows, 5a = A5. The 
fragmentation fractal dimension in nucleus-nucleus collisions is parameterized in the form; 
^AA = eo{dN/dri) + e^p. It depends on multiplicity density and nucleon fragmentation 
dimension (e^p = 0.2). Here Cq is a fitting parameter. In this approach the energy loss of 
the scattered constituent during its fragmentation in the inclusive particle is proportional 
to the value {1 — ya)- We would like to emphasize once more that dimensionless quantities, 
like momentum fractions, a scaling function and similarity parameter, are expressed via 
dimensional experimental measurable quantities. 

Figure [5] shows scaling function 'ip{z) as a function of self-similarity parameter z. The 
data z presentation corresponds to spectra of negatively charged particles produced in 
Au+Au collisions at energies yJsjsiN = 7.7 — 200 GeV in the range \ri\ < 0.5 for most 
central events (0-5%) as a function of the transverse momentum. We observe "collapse"of 
the data onto a single curve. The solid line is a fitting curve for these data. The derived 
representation shows the universality of the shape of the scaling curve for the BES 
energies. The found regularity (the shape of the function and its scaling behavior in 
the wide kinematic range) can be treated as manifestation of the self-similarity of the 
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Phc. 6: Model parameters 6 and eg as a function of collision energy y/sj^. 


structure of colliding objects, interaction mechanism of their constituents, and processes 
of fragmentation of these constituents into real registered particles. As we can see from 
Fig. [5l the scaling function exhibits two kinds of behavior; one in the low-z and the other 
in the high-z region. The first regime corresponds to saturation of the scaling function 
with the typical flattening out. The second one corresponds to the power behavior of the 
scaling function pointing on self-similarity in constituent interactions at small scales. 


C. Model parameters 

The model parameters - nucleus fractal dimension 6 a, fragmentation fractal dimension 
saa, and "heat capacity"c, are determined from the requirement of scaling behavior of 
'ip{z) (see Fig. |5]) as a function of self-similarity parameter z. The fractal dimension of the 
nucleus and the fragmentation dimension are equal to 6 a = A6 and caa = ^o{dN/dri)+epp, 
respectively. Figure [6] shows the dependence of these parameters on the collision energy. 
It was found that 6 is independent of the energy for ^/SNN > 20 GeV, cq increases with 
the energy and c is independent of the energy over the range ^/snn = 7.7 — 200 GeV. 
Note that, while 6 decreases with decreased energy for ^/snn < 20 GeV, the function 
'ip{z) preserves the shape in this range. 


D. Constituent energy loss 


The nuclear modihcation factor Raa measured at RHIG at yGViv = 62.4,130 and 
200 GeV strongly shows a suppression of the charged hadron spectra at > 4 GeV/c. 
This suppression j^-|^ was one of the first indications of a strong hnal-state modihcation 
of particle production in Au+Au collisions that is now generally ascribed to energy loss of 
the fragmenting parton in the hot and dense medium. At the BES-I energies the nuclear 


modihcation factor drastically changes with a collision energy [2^. It increases with 














b) 


Phc. 7: Momentum fraction Ua for negatively charged particle production in most central (0 — 5%) 
and peripheral (40 — 60%) Au+Au collisions at mid-rapidity (|r/| < 0.5) as a function of the 
energy y/SNjsf and hadron transverse momentum pT- 


transverse momentum as the energy decreases. The study of the evolution of the energy 
loss with collision energy has relevance to the evolution of created nuclear matter, and 
can be useful for searching for signature of phase transition and a critical point. The 
measured spectra (see Fig. [3K|) allow us to estimate constituent energy loss in charged 
hadron production in Au+Au collisions at BES energies. The estimations are based on 


a microscopic scenario of particle production proposed in 28l-l30|. The energy loss of the 


scattered constituent during its fragmentation in the inclusive particle is proportional to 
the value (1 — ya). 

Figure [7] shows the dependence of the fraction ya for central (0 — 5%) and peripheral 
(40 — 60%) Au+Au collisions on the transverse momentum over a range y/s^N = 7.7 — 
200 GeV. The behavior of ya demonstrates a monotonic growth with pt- It means that 
the energy loss associated with the production of a high-p^^ hadron is smaller than that 
with lower transverse momenta. The decrease of ya with collision centrality represents 
larger energy loss in the central collisions as compared with peripheral interactions. The 
energy dissipation grows as the collision energy increases. It is estimated to be about 20% 
at ^/sTTn = 7.7 GeV and 85% at ^/s7fN = 200 GeV for ~ 4 GeV/c for central events 
(0 — 5%), respectively. These estimations are in agreement with the energy dependence of 
the model parameters (see Fig. [6]). Discontinuity of these parameters was assumed to be 
a signature of phase transition and a critical point. 











IV. SUMMARY 


In summary, we have presented the hrst STAR results for negatively charged hadron 
production in Au+Au collisions at ^/SNN = 7.7 — 200 GeV from the Beam Energy Scan 
program at RHIC. The spectra of particle yields at mid-rapidity were measured over a 
wide range of transverse momentum. The spectra demonstrate a strong dependence on 
the collision energy and centrality which enhances with pt- An indication of self-similarity 
of negatively charged particle production in Au+Au collisions was found. It is seen in the 
presentation of spectra by scaling function -0 in dependence on self-similarity parameter 
2 ; for all energies. 

The obtained spectra of negatively charged hadrons in Au+Au collisions were used to 
estimate of a constituent energy loss as a function of energy and centrality of collisions 
and transverse momentum of inclusive particle in the ^-scaling approach. It was shown 
that the energy loss increases with the collision energy and centrality, and decreases with 


Pt- 

The energy dependence of the model parameters - the fractal and fragmentation 
dimensions and "specihc heat was studied. All parameters demonstrate smooth behavior 
with energy for most central events. It is expected that the second phase of BES will allow 
us to clarify the existence of a discontinuity and strong correlation of these parameters 
with much higher statistics. 

The presented results provide an additional motivation for the Beam Energy Scan- 
II at the STAR experiment (h3 . [^ . The large and uniform acceptance and extended 
particle identihcation is suitable for a critical point search. The BES Phase-II program 
will allow for high-statistics measurements, with an extended kinematic range in rapidity 
and transverse momentum, using sensitive observables, to reveal the structure of the QCD 
phase diagram 20| . 
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